Unitary relation between a harmonic oscillator of time-dependent frequency and a 
simple harmonic oscillator with and without an inverse-square potential 
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The unitary operator which transforms a harmonic oscillator system of time-dependent frequency 
into that of a simple harmonic oscillator of different time-scale is found, with and without an inverse- 
square potential. It is shown that for both cases, this operator can be used in finding complete sets 
of wave functions of a generalized harmonic oscillator system from the well-known sets of the simple 
harmonic oscillator. Exact invariants of the time-dependent systems can also be obtained from the 
constant Hamiltonians of unit mass and frequency by making use of this unitary transformation. The 
geometric phases for the wave functions of a generalized harmonic oscillator with an inverse-square 
potential are given. 
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It is certainly of importance to find a complete set of 
wave functions for a system of the time-dependent Hamil- 
tonian [1-17]. It has long been known that a harmonic 
oscillator of time-dependent frequency with or without 
an inverse-square potential is the system of practical ap- 
plications (see e.g. Ref. ||), where the wave functions 
are described in terms of solutions of classical equation 
of motion of the oscillator without the inverse-square po- 
tential [§-|| ■ In Ref. H it has been shown that, for a 
generalized harmonic oscillator system, the kernel of the 
system is determined by the classical action. This is one 
of the basic reasons of the fact that the wave functions 
are described by the classical solutions. On the other 



hand, it has long been noticed that there exist classical 
canonical transformations which relate the (driven) har- 
monic oscillators of different parameters (see e.g. Refs. 
|,[0|). Recently, in Ref. (ll), it has been shown that a 
driven harmonic oscillator of time-dependent frequency 
is related, through canonical transformations, to the sim- 
ple harmonic oscillator of unit mass and unit frequency 



but with a different time-scale 12 1Q]. This fact has been 
used to find the wave functions of a driven system which 
exactly agree with the known results |],[r|. 

In this Rapid Communication, we will show that for 
both oscillators with and without an inverse-square po- 
tential, there is a unitary operator which transforms the 
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harmonic oscillator systems of time-dependent frequency 
into those of the unit-mass and unit-frequency oscilla- 
tors with different time-scales. This unitary operator can 
be used to find complete sets of wave functions of the 
systems with time-dependent parameters from the well- 
known sets of wave functions of the simple harmonic os- 
cillator with [|l4j or without an inverse-square potential. 
It has been known that Q |s) there exist exact invariants 
in the systems of time-dependent parameters which have 
long been used to find the wave functions Ji] ||,[l5] |l7| . 
As might have been implied by the classical treatments 
through canonical transformations, it will also be shown 
that, the exact quantum-mechanical invariants in oscil- 
lator systems of time- dependent parameters can be ob- 
tained from the constant Hamiltonians of unit mass and 
frequency (which, certainly, are invariants in their sys- 
tems respectively), through the unitary transformation 
given here and those in Refs. ]9|,|T3|]. 

The unit-mass harmonic oscillator of time-dependent 
frequency wo(t) is described by the Hamiltonian 



P 2 , w o(t) 2 



H (x,p,t) = — , } 



with the classical equation of motion 



(1) 



x + w^(t)x = 0. 



(2) 



If we denote the two linearly independent solutions of 
Eq. (2) as Uo(t) and Vo(t), the po(t) defined by po{t) = 



\J Uq + v 2 should satisfy 
d 2 



dt 2 



SI 2 

Po + wl(t)p j 

Pa 



0. 



(3) 



with a time-constant Qq (= vqUq — uqVo). Without losing 
generality we assume that fio is positive. The wave func- 
tions ipn( x > t) 01 the system should satisfy the Schrodinger 
equation 



O (t)i/>°{x,t)=0, 



(4) 



where Oo(t) = —ih-^ + Ho(x,p,t). For the simple har- 
monic oscillator system of unit mass and frequency whose 
time is r, the wave functions tj)^{x,T) should satisfy 



O s (r)r n (x,r)=0, 



(5) 



where O s (t) = -ih£ + H s with H s = \{p 2 + x 2 ). If t 
and r is related through the relation 



dr = —^-dt, 
Po 



by defining the unitary operator 



(6) 



U w0 {p ,n ) = exp(^-x 2 )exp[--^-(ln-j^=)(xp + px)), 



(7) 



one may find the relation 



An 
Pi 



U wQ O s (T)Ul Q \ T =r(t)= O (t). 



(8) 



In Eq. (7), the overdot denotes the differentiation with 
respect to time t, while in Eq. (8) the notation " | T=T (j) " 
is to mention that r should be replaced by the function 
of t satisfying the relation (6). In a different vein, the 
relation (8) has also been noticed in Ref. p8[ . Eqs. (5,8) 
imply the following relation in wave functions; 



ip°(x,t) = U w0 i) s n \ T=T{t) . 



(9) 



As is well-known Q, the simplest choice of {^n\ n 
0, 1, 2, • • ■} may be given as 

1 



C(^ T )|r=T(t) 



\/2™n!\ 


/ irh 




X 


x exp[— 




1 




\/2 n nK 


/ irh 




x 2 


x exp[— 


2h 



fn("^ x )lr=r(t) 

U (t) - iv (t) +1/2 
Po(t) ' 

- ic ]H n (-j=x), 



(10) 



(11) 



where Co is an arbitrary real number which will be set to 
zero from now on. In obtaining Eq. (11), we make use 
of the fact: 



, fi u a -iv p 

ar = —^dt = i\ jar. 

Po u o - ivq po 



(12) 



In order to find a general expression of ipn( x i 0> we con ~ 
sider another unitary transformation. By defining S Ul (t) 
through the relations 



c 1 2 2 1-2 



(13) 



where u\ is a linear combination of Uo(t) and wo(r)> onc 
may find that the unitary operator U f given as ||l3|| 



% % 
U f = exp[-(iiix + 5 Ul (t)]exp(--uip) (14) 



satisfies the following relation 

U f O Uj = O a . 



(15) 



Therefore, the wave functions i/'n satisfying Schrodingcr 
equation of Eq. (4) may in general be written as 

1p°(x,t) = UfU w0 tfj^(x,T) | T=T ( t ) 

= 1 /^0n1/4/ "0(0 -^0(0 ^+1/2 

y/2"n\pa(t) W poW 
x exp[-(ui(r)x + 5 Ul (r))] 



*— )] 



2a v p 2 (r) po 



(16) 



a p (i) 

This wave function, of course, agrees with the known one 
P, fU^]Tl| if we consider u± as a (fictitious) particular so- 
lution. If ui = 0, the wave function given in Eq. (16) 
also agrees with that in Refs. |7|Jl5| |l7t ■ 

It may be interesting to find that how many free pa- 
rameters are in the wave function ^(x, t). First, there 
are two parameters in determining ui(t). In the case of 
til = 0, one may think that there are four parameters 
which come from determining uo(t),Vo(t). However, one 
of them is not a free parameter, since the wave func- 
tions are invariant under the multiplication of uo(t) and 
Vo(t) with same constant factor. For the simple harmonic 
oscillator of time-translational symmetry, one of the re- 
maining three parameters of u\ = is simply related to 
a time-shifting of the wave function. This can be seen 
from the fact that, for the unit frequency case, the uq 
and Vq can be taken as cos(i + to) and C sin(r + [3 + to), 
respectively, with real constants C, (3, ta- 
li one considers p s (t) satisfying 

d 2 n 2 



-Ps + p s 

dr 2 pi 



0, 



(17) 



and a simple harmonic oscillator of unit mass and fre- 
quency and with time r' which is related to r as 



dr' = —^-dr, 

A 



by defining 
U s = exp( 



dT ^) eM -^ { ln-E^){xp + px)), (18) 



1 

2hp s ~ 2~h" 



one may find that 



Pi 



U s O s (r')Uj | t , =t , (t) =O s (t). 



(19) 



The wave functions Vn( T ) defined by 

4>n( T ) = U s i> s n (x,T') | r / =T / (T ) 

then satisfy the Schrodinger equation O s (r)'0^ = 0; In 
fact, ipn(r) is closely related to the wave functions of the 
squeezed states p9|j20| ,P|. 

One may think that a more general expression of the 
unitary operator, U w o, may be obtained by combining 
use of U w o and U s . This, however, is not the case as can 
be seen from the relation 

U w0 (po, O, )U s | r=T ( t )= U w0 (poPs, ^ £l s ) \r=r(t), (20) 

which is in accordance with the number counting of free 
parameters in ip^(x,t). 

The harmonic oscillator of unit mass and frequency 
with an inverse-square potential is described by the 
Hamiltonian l]14| 



Hf. 



2 2 

p z x g 



(21) 



We only consider the case of g > — Ti /8, and the region 
of x > 0. By defining a = ±(1 + 8g/h 2 )^ 2 and 



or 



(22) 



the wave function (j) s n satisfying the Schrodinger equation 
Of{T)<t> s n = is given as @ 

_ / 4 -,i/4/ r(n+l) 1/2 i(a n+a+ i) r 
V y T(n + a + V 



1 11 

x(^)(-«)/«exp(-^)L»(^,. 



(23) 



By defining O™ as 



9 p 



oT(t) = -mf t + | + »„ 2 (t)y + ±, (24) 

as in the case without the inverse-square term, one may 
find the relation 



^U w0 OT(r)Ul \ T=T(t) = 0™(t). (25) 
Po 



In deriving Eq. (25), we make use of the commutator 
relation 



[xp + px, —] = 4ih—. 



(26) 



For a further generalization, we define a unitary operator 



CT fl = exp[-(Max 2 - — )a; 2 ] exp[i^-(xp + px)}, (27) 



where M is a positive function of i, and a(t) is a real 
function. One may then easily find the relation 



U g O™Ui = O l 



at 



(28) 
(29) 



where (see Ref. |J) 



Hm = 2M(t) ~ a ^ x 'P + P x 
+±M(t)c(t)x 2 



9 1 

M(i) x 2 



(30) 



with 



s 1 <fy/M , Id, 

c W ^) + ^— + 4a 2 -2--(Ma). 



For convenience fl3|| , we consider the equation 



— (Mx)+w 2 (t)x = 0, 



(31) 



where w 2 (t) = w%(t) + -±=^-0-. The two linearly in- 
dependent solutions u[t),v(t) of Eq. (31) can be given 
from Uo(t),Vo[t) as u(t) — -^2=,v(t) = -^=, so that one 
may find the relation Slo = M(vu — iiv). We also define 
the pit) as pit) — -^=. The wave function </>„ of the sys- 
tem described by the Hamiltonian Hi n {x 1 p 1 1) can then 
be obtained as 



(32) 



<p n = U G (j) S n iT) | T=T ( t ) 

l Rp2 ; T(n + a + l j 
^ u-w . (2w+a+1 ^ O a; w 2 a + iv4 
p ftp 2 

X 6Xph S ( ^ ~ <M £ ~ 2iM «)]^(^)' (33) 



where 



Ug = U„U w o- 



(34) 



For a = 0, the wave functions <pn agree with those in 
Refs. 0,0- As in Ref. ||], by considering the kernel of 
the system , it may be easy to see that the wave func- 
tions <j>n{x, t) form a complete set. The form of <p n in Eq. 
(33) indicates that, even for the system described by the 
constant Hamiltonian H? n given in Eq. (21), there are 
wave functions whose probability density distributions 
pulsate as in those of the squeezed states. 

For the system of the Hamiltonian i?j n , if M(t), u^j(t), 
and a(t) are periodic with a period T, one may study the 



non-adiabatic geometric phases (21 22 1. The wave func- 
tion 4> n is (quasi)periodic, only if pit) is periodic. The 
condition for periodic pit) with the period T' (= T or 
2T) has been analyzed in Ref. {23l. Here, we only con- 



sider the case of such a periodic pit). The overall phase 
change of <f> n under the T' evolution is given as 



Xn = -(2n + a + l) 



T' 



n 



dt. 



la m) 

The expectation value of the Hi n can be evaluated by 
making use of the relation 



d d 
H in <t>n = (in—U G )(f> s n + U G ih—<, 



(35) 



From the fact that 



^ = i 4^ = ^ + a + 1 ^^ (36) 
one may find that the geometric phase 7„ for the wave 
function <f) n under the T' evolution is written as 
1 f T ' 

ln=Xn + T {4>n\H in \(j) n )dt 

= 5?T / ^ M P 2 + 2Mapp)dt(cj> s n \x 2 \<f> s n ) 

= (2n + a + l)-j- / iMp 2 + 2Mapp)dt. (37) 
"o Jo 

The unitary operators can be used in finding the exact 
invariants for the cases without and with the inverse- 
square potential from H s and Hf n , respectively. First 
of all, it is clear that H s and Hf n are invariants in the 
systems they describe, respectively. For the system de- 
scribed by Hoix,p,t), if we only consider the case of 
u\ = 0, the invariant Jo is obtained by applying the uni- 
tary transformation to the invariant H s 



Iq — U w qH s U wQ 

= JrK— ) 2 + (pop-px) 2 ) 

21 Iq po 



(38) 



which agrees with those in Refs. @-@,[l5|,[l6| . For the sys- 
tem described by H in ix,p,t), the invariant is again given 
from the invariant Hf„ as 



I m = U G Ht n Ul 
2fi P 

+ {pp - (Mp + 2Map)x} 2 + 2p 2 ^]. (39) 

For the case of a = 0, the invariant ij n reduces to the 
known one Jl7|] . One can explicitly check that the invari- 
ant Ii n indeed satisfies the relation 

dl 

in^ + [I in ,H in }=0. (40) 

Alternatively, making uses of Eqs. (35,36) and relying on 
the completeness of the set {</>^|n = 0,l,2,---},a simple 
proof of Eq. (40) may also be possible. 

In summary, we have found a unitary operator 
which transforms a harmonic oscillator system of time- 
dependent frequency into that of a simple harmonic os- 
cillator of different time-scale, with and without the 
inverse-square potential. Making use of the unitary op- 
erator, the exact invariants and wave functions of the 
time-dependent systems have been evaluated from the 
well-known results in the corresponding system of con- 
stant Hamiltonians. It should be mentioned, however, 
that the classical solutions of the time-dependent har- 
monic oscillator system must be found for actual appli- 
cations, while the classical equation (see Eq. (2)) is for- 
mally equivalent to a one-dimensional time-independent 
Schrodinger equation (of arbitrary potential). The classi- 
cal correspondent of unitary transformation is the canon- 
ical transformation which has been studied in the model 
Jjl] , [ic| . It would be interesting if the relationship could 
be used in finding relations among the quantities in clas- 



sical and quantum mechanics such as that between the 
geometric phases and the Hannay's angle (see Ref. p3[). 
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